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Hydrogen adatoms are shown to generate magnetic moments inside single layer graphene. Spin 
transport measurements on graphene spin valves exhibit a dip in the non-local spin signal as a 
function of applied magnetic field, which is due to scattering (relaxation) of pure spin currents by 
exchange coupling to the magnetic moments. Furthermore, Hanle spin precession measurements 
indicate the presence of an exchange field generated by the magnetic moments. The entire experi- 
ment including spin transport is performed in an ultrahigh vacuum chamber, and the characteristic 
signatures of magnetic moment formation appear only after hydrogen adatoms are introduced. Lat- 
tice vacancies also demonstrate similar behavior indicating that the magnetic moment formation 
originates from p z -orbital defects. 

PACS numbers: 72.80.Vp,85.75.-d,75.30.Hx,72.25.Rb 



Many fascinating predictions have been made regard- 
ing magnetism in graphene including the formation of 
magnetic moments from dopants, defects, and edges [1- 
8]. While several experimental techniques provide in- 
sight into this problem [9-21], lack of clear evidence for 
magnetic moment formation hinders development of this 
nascent field. Studies based on bulk magnetometry [9- 
16] directly measure magnetic properties, but because it 
measures the total magnetic moment (not just the sig- 
nal from graphene) it is difficult to rule out artifacts 
from environmental magnetic impurities. Transport [17- 
19] and scanning tunneling microscopy (STM) [20, 21] 
locally probe the graphene, but so far these measure- 
ments have been charge-based, so data are subject to var- 
ious interpretations [22]. Thus, in order to convincingly 
demonstrate the formation of magnetic moments inside 
graphene due to dopants and defects, it is essential to 
employ techniques that directly probe the intrinsic spin 
degree-of-freedom of the magnetic moment while ensur- 
ing the signal originates from the graphene sheet under 
investigation. 

In this Letter, we utilize pure spin currents to demon- 
strate that hydrogen adatoms and lattice vacancies gen- 
erate magnetic moments in single layer graphene. Pure 
spin currents are injected into graphene spin valve de- 
vices and clear signatures of magnetic moment formation 
emerge in the non-local spin transport signal as hydro- 
gen adatoms or lattice vacancies are systematically in- 
troduced in an ultrahigh vacuum (UHV) environment. 
Specifically, introduction of these point defects generate 
a characteristic dip in the non-local signal as a func- 
tion of magnetic field. This feature is due to scattering 
(relaxation) of pure spin currents by localized magnetic 
moments in graphene and is explained quantitatively by 
a phenomenological theory based on spin-spin exchange 
coupling between conduction electrons and magnetic mo- 
ments. Furthermore, we observe effective exchange fields 



due to this spin-spin coupling, which arc of interest for 
novel phenomena and spintronic functionality [23-26] but 
have not been seen previously in graphene. Thus, these 
results provide the most clear and direct evidence for 
magnetic moment formation in graphene and demon- 
strate a method for utilizing localized magnetic moments 
to manipulate conduction electron spins. 

For a systematic investigation, the spin transport mea- 
surement is first performed on a pristine single layer 
graphene (SLG) spin valve as a control measurement. 
Then, dopants/defects are controllably introduced to the 
SLG and the measurement is repeated. The sample re- 
mains in UHV during the entire process. Therefore, 
observed signatures of magnetic moment formation are 
caused by the adsorbed hydrogen or lattice vacancies. 

Experiments arc performed on non-local SLG spin 
valves [27-29] (Fig. la) consisting of two outer Au/Ti 
electrodes (a and d) and two ferromagnetic (FM) Co elec- 
trodes that make contact to SLG across MgO/Ti0 2 tun- 
nel barriers (b and c) . The Co electrodes are capped with 
5 nm AI2O3 to protect from hydrogen exposure. The tun- 
nel barrier and capping layer are present only at the site 
of the FM electrodes, leaving the rest of the graphene 
uncovered. The device is fabricated on a Si02/Si sub- 
strate (300 nm thickness of Si0 2 ) where the Si is used as 
a back gate. Details of device fabrication are published 
elsewhere [29]. 

The charge and spin transport properties of pristine 
SLG spin valves are measured at 15 K using lock-in tech- 
niques. The gate dependent resistivity (pa) of a repre- 
sentative sample A (black curve in Fig. lb) exhibits a 
maximum at the gate voltage (Vg) of V, which defines 
the Dirac point (Vb = V). This sample exhibits mobil- 
ity (//) of 6105 cm 2 /Vs. To investigate spin transport in 
the SLG device (Fig. la), a current (7) is applied between 
electrodes b and a, injecting spin-polarized carriers into 
graphene directly below the FM injector, b. The spin 
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FIG. 1. The effect of hydrogen exposure on charge and spin 
transport in SLG at 15 K. (a) Schematic illustration of the 
non-local spin valve device, (b) Gate dependent resistivity 
for the pristine graphene (black) and following exposure to 
atomic hydrogen for 2 s (red) and 8 s (blue). Upon hydrogen 
doping, the Dirac point shifts from V to -1 V. (c) Non-local 
spin transport measurement for pristine graphene. (d) Hanle 
spin precession measurement on pristine graphene. (e,f ) Non- 
local spin transport measurements after atomic hydrogen ex- 
posure for 2 s and 8 s, respectively. Both curves exhibit a dip 
in Rnl at zero applied field, which is caused by spin relax- 
ation induced by localized magnetic moments. 



AR NL of 8.8 fi (sample A with V G - V D = -15 V). A con- 
stant spin-independent background is subtracted from all 
Rnl data presented in this study. Out-of-plane mag- 
netic fields are applied to generate spin precession, and 
the resulting data (Fig. Id, red for parallel, blue for an- 
tiparallel) are fit by the standard Hanle equation [28, 29] 
(solid curves) to determine the spin lifetime (t so = 479 
ps) and diffusion coefficient (D — 0.023 m 2 /s ). The cor- 
responding spin diffusion length is A = v Dt so — 3.3 pm. 
Based on these values and a non-local spin signal of 8.8 
f2, the spin polarization of the junction current (Pj) is 
calculated to be 20% [30, 31]. 

Atomic hydrogen is introduced to spin valve devices at 
15 K at a chamber pressure of 1 x 10 -6 torr [31]. Fol- 
lowing 2 s hydrogen exposure, the gate dependent p G 
(red curve in Fig. lb) is dramatically increased. An ad- 
ditional 6 s of exposure (8 s total) further increases pc 
(blue curve of Fig. lb) and decreases the mobility to 495 
cm 2 /Vs. Based on the change in the resistivity, we make 
an order of magnitude estimate for the hydrogen cover- 
age of 0.1% [31]. Accompanying the changes in charge 
transport are also changes in spin transport. Figures le 
and If display Rnl of sample A at Vq — Vd = -15 V as 
a function of B apps following 2 s and 8 s of exposure, 
respectively. The initial ARnl of 8.8 fl is reduced to 
2.6 O after 2 s of hydrogen exposure and further reduced 
to 1.4 il after 8 s. Interestingly, the Rnl scans exhibit 
a dip centered at zero applied field. The dip in Rnl is 
prevalent for both up and down sweeps of B apPiy at all 
measured gate voltages and has been reproduced on mul- 
tiple samples following hydrogen exposure. The ratio of 
the dip magnitude to ARnl is found to increase with in- 
creasing hydrogen exposure (comparing Fig. le and If), 
indicating the dip feature is dependent on the amount of 
adsorbed hydrogen. 

To understand the origin of the dip in Rnl, we exam- 
ine the expression for non-local resistance generated by 
spin transport [30], 
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population diffuses along the sample as a pure spin cur- 
rent (a;-axis) and the spin density is measured at the FM 
spin detector, c, as a voltage difference (V) between elec- 
trodes c and d. An applied magnetic field (B appy ) along 
the electrode magnetization direction (y-axis) is used to 
control the relative orientation of spin injector and detec- 
tor magnetizations. For parallel alignment, the measured 
non-local resistance (Rnl — V/I) is positive whereas for 
antiparallel alignment Rnl is negative. The non-local 
spin signal is defined as the difference between parallel 
and antiparallel states (ARnl = Rnl ~ R-m, )• A typi- 
cal scan of Rnl as a function of B app _ y (Fig. le) displays 
discrete jumps as the electrode orientation changes be- 
tween parallel and antiparallel. This sample exhibits a 
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where R G = p G ^/w is the spin resistance of graphene, 
w is the graphene width, Rp = ppXp/Aj is the spin re- 
sistance of the cobalt, pp is the cobalt resistivity, A_f is 
the cobalt spin diffusion length, Aj is the junction area, 
Pp is the spin polarization of cobalt, i?i and R2 are the 
contact resistances of the spin injector and detector, re- 
spectively, and L is the distance from injector to detector. 
This equation shows that the spin density at the detec- 
tor electrode depends on both charge and spin proper- 
ties. First, we confirm that the SLG resistivity does not 
change with magnetic field, so the dip is not related to 
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FIG. 2. (a) A minor loop scan shows that the dip in Rnl 
for parallel alignment (black) becomes a peak for antiparallel 
alignment (red), indicating the feature is due to increased spin 
relaxation, as opposed to an artifact of the background level, 
(b) Fitting the dip in Rnl based on the model of spin relax- 
ation by paramagnetic moments (data in black, fit in red), (c) 
Field dependence of longitudinal (red) and transverse (black) 
spin lifetimes, (d) Hanle precession data following 8 s hydro- 
gen exposure (red) is fit using equation 3 (black curve). 



changes in charge transport [31]. Second, we verify that 
the dip is not related to hydrogen-induced changes to 
the magnetic properties of the FM electrodes. Specifi- 
cally, the effect of hydrogen exposure is reversible upon 
thermal cycling to room temperature and the anisotropic 
magnetoresistance of the Co electrodes are not affected 
by hydrogen exposure [31]. Next, we perform minor loop 
analysis on sample B (Fig. 2a) by reversing the magnetic 
field sweep immediately after the first magnetization re- 
versal. The inversion of the dip in the antiparallel state 
(red curve) proves that the dip is due to increased spin re- 
laxation at low fields. Furthermore, we rule out hyperfine 
coupling to nuclear spins as the origin of this increased 
spin relaxation [31]. 

As we discuss in the following, emergence of the dip 
following hydrogen adsorption identifies magnetic mo- 
ment formation in graphene. The dip in Rnl is a char- 
acteristic feature of spin relaxation from exchange cou- 
pling with localized magnetic moments, and can be il- 
lustrated from a simple textbook example of two cou- 
pled spins in a magnetic field. The Hamiltonian is given 
by H = A ex S e ■ Sm + g e ^BS e • B app + qm^bSm • B app , 
where S e is the conduction electron spin, Sm is the spin 
of the magnetic moment, g e and qm are the respec- 
tive ^-factors, and A ex is the exchange coupling strength 
[32, 33]. Due to the presence of the exchange coupling, 



the individual spins are not conserved; only the total spin 
Stot — S e + Sm is conserved. For the case where both 
S e and Sm are spin- A, the quantum mechanical eigen- 
states in zero magnetic field are the well-known singlet 
{Stot — 0) and triplet {S to t = 1) spin states [34]. At 
higher magnetic fields the Zeeman terms dominate and 
the two spins decouple so that the magnitudes and z- 
components of S e and Sm become good quantum num- 
bers, similar to the Paschen-Back effect [34]. Thus, 
the dip in Rnl is qualitatively explained by the non- 
conservation of at low fields due to the presence of 
exchange coupling with magnetic moments. 

To quantitatively analyze the experimental data, we 
must consider that a conduction electron will interact 
with many localized magnetic moments. Thus, the terms 
in the Hamiltonian involving the conduction electron are 
given by H e = r]MA ex S e -{S M )+ge^BSe-B app = g e ^BS e - 

{ji ex + Bapp^j where t)m is the filling density of magnetic 

moments. The averaging (• • •) is over the ensemble of 
magnetic moments and the effective field generated by 
the exchange interaction is B ex = wA^Lj^l As the 
spins diffuse through the lattice they experience varying 
magnetic moments which results in varying Larmor fre- 
quencies. In the local frame associated with the electrons 
this can be described by a time-dependent, randomly 

fluctuating magnetic field, B ex {t) = B ex + AB ex (t). For 
the Rnl measurements, the longitudinal spin relaxation 
due to a fluctuating field is given by [35] , 
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where AB is the rms fluctuation and r c is the corre- 
lation time [31]. The spin relaxation rate due to the 
exchange field is described by a Lorentzian curve which 
depends explicitly on the applied field, B 
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result- 
ing in strong spin relaxation at low fields and suppressed 
spin relaxation at high fields. Due to the presence of 
B ex ,y in equation 2, ferromagnetic ordering will produce 
a dip in Rnl that is centered away from zero and is 
hysteretic, while paramagnetic ordering will produce a 
non-hysteretic dip centered at zero field. Thus, the mag- 
netic moments measured in these experiments are para- 
magnetic. The total longitudinal spin lifetime, Tl° tal , 
of conduction electrons is dependent on both the usual 
spin relaxation due to spin orbit coupling (t so ) and lon- 
gitudinal spin relaxation from the exchange field (T-f x ), 
such that {T{ otal )~ 1 = (rf*) -1 + (r so ) _1 . We apply the 
above model to the non-local spin transport data pre- 
sented in Fig. If (sample A) and fit using equation (1), 
A = BT{ otal , and equation (2) [31]. The resulting fit 
(red line in Fig. 2b) replicates the shape and magnitude 
of the dip measured in Rnl (black line in Fig. 2b). The 
field dependent T[ otal (Fig. 2c), exhibits a minimum of 
464 ps at zero field and increases asymptotically towards 



4 



531 ps for large B a 



The values obtained for 



AB and r c are 6.78 mT and 192 ps, respectively. The 
field-dependent spin relaxation following atomic hydro- 
gen exposure, which emerges as a dip in Rnl, is a clear 
signature of paramagnetic moment formation. 

Spin precession measurements provide further evidence 
for the presence of magnetic moments. Figure 2d shows 
spin precession data for sample A (8 s exposure, Vq — Vd 
= -15 V) with FM electrodes in the parallel alignment 
state. The Hanle curve has considerably narrowed com- 
pared to the precession measurements obtained prior to 
hydrogen adsorption (Fig. Id). The sharpening of the 
Hanle curve results from the presence of an exchange 
field. The injected spins precess around a total field 
Btot = B apPyZ + B ex ^ z (along z-axis) that includes not 
only the applied field, but also the exchange field from 
the paramagnetic moments. At 15 K and B appz < 100 
mT, the magnetization is proportional to the applied 
field so that B ex ^ z = kB appz , where k is a proportion- 
ality constant. Thus, the spins precess about B to t with 
frequency w = g e \i B B to tl^ = 5e(l + k)n B B app , z /h = 
gln B B app ^ z /Ti. To properly account for the enhanced 
g-factor induced by the magnetic moments, the Hanle 
equation must be modified to 
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is the transverse spin lifetime. As shown in 
Fig. 2c, the T\ otal is related to, but different from T[ otal 
[31] . Using the field dependent T*°* oZ , the precession data 
(red circles of Fig. 2d) is fit to equation 3 (black line) to 
yield a value of g* e = 7.13. Physically, g* e > 2 corre- 
sponds to an enhanced spin precession frequency result- 
ing from the exchange field. A detailed discussion of the 
Hanle fitting and the gate-dependent properties of the 
exchange field are provided in the supplementary infor- 
mation [31]. The dramatic narrowing of the Hanle peak 
combined with the emergence of a dip in Rnl provides 
the most direct evidence to date for the formation of 
magnetic moments in graphene due to the adsorption of 
atomic hydrogen. 

We now turn our attention to lattice vacancy defects in 
graphene. Several theoretical works suggest the similar- 
ity of magnetism due to vacancies and hydrogen-doping 
[1, 4], as both should create magnetic moments inside 
graphene due to the removal/hybridization of p z orbitals. 
To produce lattice vacancies in pristine SLG spin valves, 
we perform Ar-sputtering at low energies and examine 
the subsequent non-local spin transport. We again ob- 
serve the emergence of a dip in Rnl and narrowed Hanle 
curve, indicating the formation of paramagnetic moments 
in graphene [31]. Given the very different chemical and 
structural properties of lattice vacancies compared to ad- 
sorbed hydrogen, the observation of similar features in 
the spin transport data provide strong evidence that the 



magnetic moments are created by the removal of p z or- 
bitals from the 7r-band, as predicted theoretically. 

In conclusion, clear signatures of magnetic moment for- 
mation are observed in both the non-local spin transport 
and Hanle precession data, which emerge only after expo- 
sure to atomic hydrogen or lattice vacancies. The results 
and techniques presented here are important for future 
developments in magnetism and spintronics. 
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I. ATOMIC HYDROGEN EXPOSURE 

A commercial Omicron source is used to expose SLG 
spin valve devices to atomic hydrogen at 15 K. Diatomic 
hydrogen is cracked inside a tungsten capillary tube that 
is heated by electron bombardment. The amount of hy- 
drogen introduced to the chamber is controlled via a leak 
valve, which is tuned to maintain a chamber pressure of 
1 x 10~ 6 torr (the base pressure of the chamber is below 
1 x 10 -9 torr). The heating power of the Omicron source 
is determined by the high voltage (HV) applied to the 
capillary and the emission current between capillary and 
filament (I cm ). We use the parameters HV=1 kV and 
Iem=80 mA. The distance from source to sample is 100 
mm. A shutter positioned between the SLG spin valve 
and hydrogen source is used in order to control the expo- 
sure time. Additionally, deflector plates are used to steer 
any charged ions away from the sample. 

As discussed in the main text, the exposure of SLG 
spin valves to atomic hydrogen substantially modifies 
charge transport properties, such as pa and p (Fig. 1 of 
main text). We obtain an order of magnitude estimate 



for the hydrogen concentration based upon the changes 
in charge transport properties assuming adsorbed hydro- 
gen induces resonant scattering. Comparing with previ- 
ous experimental work on resonant scattering in graphene 
via fluorine doping 1 and lattice vacancies 2 , the hydrogen 
concentration is estimated to be on the order of 0.1% 
for 8 s hydrogen exposure to sample A. This indicates 
samples are in the dilute limit of hydrogen coverage. 



II. CONTROL EXPERIMENTS TO 
DETERMINE THE ORIGIN OF THE DIP IN R NL 

A. Excluding field-dependent resistivity effects 

Equation (1) of the main text shows that changes in 
the resistivity of graphene could change the spin density 
at the detector electrode. Thus, magnetic field depen- 
dent changes to the graphene resistivity could in princi- 
ple produce variations in Rnl- To investigate this pos- 
sibility in hydrogen-doped graphene, gate dependent pc 
measurements are performed on sample B at three dis- 
tinct in-plane magnetic fields (Fig. 1). The three curves 
are indistinguishable, showing that the applied magnetic 
field has no effect on the measured pa of hydrogen-doped 
samples. Therefore, the dip in Rnl does not originate 
from field-dependent changes in resistivity. 



B. Excluding nuclear spin effects 

We investigate the possibility that the observed dip 
in the non-local resistance originates from hyperfine cou- 
pling between conduction electrons and nuclear spins. In 
graphene, this scenario is unlikely due to the small abun- 
dance of intrinsic nuclear spins in carbon (>98% of car- 
bon is 12 C, which has no nuclear spin) and a lack of con- 
tact hyperfine coupling in the p 2 -orbitals that make up 
the conduction and valence bands. Nevertheless, investi- 
gating the situation is necessary because the adsorption 
of hydrogen on graphene may alter the hyperfine cou- 
pling. Two effects that could in principle alter the non- 
local resistance include hyperfine coupling to dynamically 
polarized nuclear spins 3 ' 4 and organic magnetoresistance 
(OMAR) 5 - 6 . 

The effect of OMAR has previously been observed in 
carbon Ceo and functionalized carbon-based polymers 5,6 . 
OMAR originates from hyperfine induced spin mixing 
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FIG. 1. Investigating the source of the dip measured in Rnl- 
Gate-dependent pa of SLG after hydrogen exposure measured 
at several discrete in-plane applied magnetic fields. The ap- 
plied field has no effect on the charge transport behavior of 
the graphene device. 
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FIG. 2. Test of nuclear spin effects. The non-local spin 
transport signal is measured for 10 min while -B a pp,j/=30 mT 
(black). This ensures any possible nuclear spins will reach 
steady state. The field is ramped to zero, and held constant 
for the following 20 min (red). The observed response ex- 
cludes hyperfine coupling to nuclear spins as the origin of the 
dip in Rnl- 



between singlet and triplet states and manifests itself as a 
magnetic field dependent resistivity. However, as shown 
in Fig. 1, resistivity does not change as a function of 
applied in-plane magnetic field. This confirms OMAR is 
not responsible for the observed dip in Rnl- 

The effect of dynamic nuclear polarization (DNP) was 
demonstrated clearly by Salis et al. 4 and Chan et al. 3 
who investigated GaAs non-local spin valves at low tem- 
peratures. Specifically, Salis et al. observed a dip in Rnl 
at zero applied magnetic field, similar to the dip we ob- 
serve in hydrogen-doped graphene. They attributed their 
dip to precessional spin dephasing caused by hyperfine 
coupling to dynamically polarized nuclear spins. To de- 
termine whether such nuclear spin effects are present in 
hydrogen-doped graphene we perform a series of tests. 

Test 1: Nuclear spin relaxation times are typically 
long (^minutes), and therefore, slow dynamics at this 
time scale are a characteristic of effects related to DNP. 



This is manifested in non-local spin transport data as 
a "lab time" dependence 4 . fn our investigation of 
hydrogen-doped graphene, we do not observe a lab time 
dependence or a magnetic field ramp rate dependence. 

Test 2: A characteristic feature of hyperfine coupling 
to the nuclear spin bath through DNP is the nuclear 
field's linear dependence on the applied field 4 . Specifi- 
cally, at zero field the nuclear spin bath depolarizes slowly 
over time. The depolarization of the nuclear spin bath is 
evident in Rnl data as the gradual decrease and even- 
tual disappearance (after a few minutes) of the dip when 
the applied field is set to zero (Figure 1 of Salis et al. 4 ). 
We perform this test on hydrogen-doped graphene spin 
valves, as shown in Fig. 2. First, the magnetic field is 
held at -30 mT for 10 minutes to ensure any possible 
nuclear spin transients reach steady state. The Rnl is 
measured continuously during this period and exhibits 
a value of ^0.92 ft. Then the magnetic field is quickly 
reduced to zero, coinciding with the immediate drop of 
Rnl to ^0.80 fi. This drop occurs because B app y =0 
mT is at the center of the dip in Rnl- The Rnl is 
measured over the following 20 minutes, for which the 
observed value remains unchanged at ^0.80 f2. This in- 
dicates that the magnitude of the dip is independent of 
lab time. If the dip were due to hyperfine coupling to 
dynamically polarized nuclear spins, then the magnitude 
of the dip would gradually decrease to zero (i.e. Rnl 
would increase back to ~0.92 il). Because this behavior 
is not observed, the dip in Rnl cannot be due to DNP. 

Test 3: For the case of DNP, applying a constant out- 
of-plane magnetic field during a non-local spin transport 
measurement (in-plane field scan) results in a charac- 
teristic feature of nuclear depolarization and repolariza- 
tion as the in-plane field crosses zero (Figure 1 of Chan 
et al. 3 ). We perform this measurement on hydrogen- 
doped graphene and observe no evidence of depolariza- 
tion/repolarization features. Together, these three tests 
show that the dip in Rnl is not due to hyperfine coupling 
to dynamically polarized nuclear spins. 

The above investigations of OMAR and DNP conclu- 
sively exclude the possibility of hyperfine coupling with 
nuclear spins as the source of the observed dip in the 
non-local spin transport. 



Excluding changes in magnetization of the 
electrodes 



Here we examine the possibility that the observed be- 
havior in Rn l after hydrogen exposure is caused by spu- 
rious changes to the magnetic electrodes. It is known 
that hydrogen adsorption onto ferromagnetic materials 
can alter their magnetic properties 7-11 . In the present 
experiment, following the standard fabrication procedure 
for tunneling contacts to graphene as discussed in previ- 
ous work 9 , the magnetic electrodes are capped with 5 nm 
AI2O3. Also, we note that in the Rnl data presented in 
this work, the coercive fields of the electrodes remain un- 
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FIG. 3. Reversibility of effects due to hydrogen doping of SLG. a, Gate dependent resistivity for pristine SLG at T=15 K. 
Data set corresponds to sample A. b, Gate dependent resistivity after 8 s hydrogen exposure at T=15 K. c, Gate dependent 
resistivity at T=15 K, measured after hydrogen exposure at 15 K followed by warming to room temperature and subsequent 
cooling to 15 K. d, Rnl for pristine SLG. e, Rnl after 8 s hydrogen exposure, f, Rnl after temperature cycling, g, Hanle spin 
precession data for pristine SLG. Blue (red) data is for parallel (antiparallel) magnetizations, h, Hanle spin precession data 
after 8 s hydrogen exposure, i, Hanle spin precession data after temperature cycling. All -Rjvl and spin precession data are 
measured at Vg — Vb=-15 V. 



changed after hydrogen exposure. 

Further it is known that the bonding energy for hydro- 
gen adsorption on transition metals is several eV (~2.6 
eV for Co 11 ), suggesting that H-Co chemisorption would 
be robust to temperature cycling. Here we note that 
the observed behavior in Rnl and spin precession data 
after the introduction of atomic hydrogen at T=15 K 
are reversible upon temperature cycling. At cryogenic 
temperatures, exposure of SLG graphene spin valves to 
hydrogen dramatically alters the charge and spin prop- 
erties. The gate dependent resistivity increases, whereas 
the magnitude of Rnl decreases and exhibits a dip at low 
field. The observed dip is due to an increase in the spin 
scattering. Also, precession measurements exhibit a nar- 
rowing of the Hanle curve. Next, the sample is warmed to 
room temperature and subsequently re-cooled to T=15 
K. The spin and charge transport properties are then 
re-measured. We have found that the charge transport 



properties of the gate dependent resistivity and mobil- 
ity nearly recover to the pristine values. For instance, 
for sample A, the mobility is /i=6105 cm 2 /Vs for pris- 
tine graphene. After 8 s hydrogen exposure the mobility 
decreases to /i=495 cm 2 /Vs. After temperature cycling, 
the mobility recovers to /i=5450 cm 2 /Vs. The recovery 
indicates the effect of hydrogen is removed through des- 
orption, cluster formation, or a combination of the two. 
After temperature cycling, the spin signal Rnl contains 
no indication of a dip near low field and has increased 
to 5.95 fl. Also, the spin precession data broadens and 
conventional Hanle analysis (equation 14) yields r so =353 
ps and £>=0.022 m 2 /s. As mentioned above, changes in 
the Co electrodes should persist after temperature cy- 
cling and so cannot explain the observed behavior. Upon 
re-hydrogenation at cryogenic temperatures the key fea- 
tures (dip in Rnl, narrowed Hanle curve) return. Fig. 3 
summarizes the changes in resistivity, Rnl, an d spin pre- 
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FIG. 4. Magnetoresistance of Co wires before and after 
hydrogen exposure at T=10 K. a, AR/R (%) of a single Co 
wire measured before hydrogen exposure (red circles) for a 
wire oriented along B ap p,y (0 °). The wires were fabricated 
using standard e-beam lithography and angle evaporation in 
the same fashion as for producing tunnel barrier electrodes on 
SLG. The observed features are due to anisotropic magnetore- 
sistance (AMR) and reflect the magnetization of the Co wire. 
No effect is observed after 20 s of hydrogen exposure (black 
diamonds), b, AR/R (%) of a Co wire oriented 45 ° with 
respect to B app ^y before hydrogen exposure (red circles) and 
after 20 s exposure (black diamonds), c, AR/R (%) of a Co 
wire oriented 90 ° with respect to B ap p,y before hydrogen ex- 
posure (red circles) and after 20 s exposure (black diamonds). 
In all cases the wires are oriented in-plane. 



cession upon temperature cycling. 

To comprehensively verify that changes in magnetism 
of the Co electrodes is not the cause of the observed dip, 
we fabricate Co wires in the same fashion as for tun- 
nel barrier contacts to graphene. Three 300 nm wide Co 
wires with equal length (L=200 /im) and 80 nm thick are 
oriented at ° , 45 ° , and 90 ° with respect to an in-plane 
applied field. Figure 4 shows the magnetoresistance plot- 



ted as AR/R (%) vs. B appy for the three different ori- 
entations before and after 20 s hydrogen exposure. We 
note that 20 s hydrogen exposure is significantly more 
than typically needed to generate observable effects in 
Rnl- As can be seen in Fig. 4, no effect is seen in the 
magnetoresistance of the Co wires upon hydrogen expo- 
sure. The switching fields and shape remain unchanged 
along with the anisotropy behavior. Therefore, we rule 
out effects of alterations in the Co wire magnetization 
as a possible source for the origin of the observed dip at 
zero field in Rnl- 



III. EFFECTIVE EXCHANGE FIELD MODEL 

To quantitatively model the experiment, we consider 
electron spins S e moving in an effective magnetic field 
of randomly positioned local magnetic moments Sm of 
filling density r)M- Each electron feels the average spin 
interaction 

9e^BS e ' B app 
= gefJ-BSe ■ (B ex + B app ) (1) 



and 



B ex — VMA ex (SM)/g e fiB 



(2) 



where B ex is the effective exchange magnetic field, g e is 
the electron g- factor, is the Bohr magneton, and A ex 
is the strength of the exchange coupling between S e and 
Sm- The averaging (...) is over the ensemble of magnetic 

moments. The effective exchange field B ex contributes to 
the Larmor frequency and enhances the electron g-factor. 

As the spins diffuse through the lattice they experience 
varying magnetic moments which results in varying Lar- 
mor frequencies. In the local frame associated with the 
electrons this can be described by a time-dependent, ran- 
domly fluctuating magnetic field, B ex (t) — B ex +AB ex (t) 
with the rms value given by the time average 



{AB™f a = {[AB ex , a {t)f 



(3) 



(AS 



rms \ 2 _ 
ex 



=(AB:r)i + (AB:r) 2 y + (AB, 



(4) 



where a is an xyz component index. The time scale of 
the fluctuation is given by a correlation time r c defined 
by 



(AB ex (t) ■ AB ex {t-t')) t ex exp(-i'/r c ) 



(5) 



Spin relaxation resulting from a randomly fluctuating 
magnetic field has been solved in the review article by 
Fabian et al. 12 , Section IV. B. 2 and is mathematically 
analogous to the D'yakanov-Perel model 13 . For the non- 
local spin signal geometry, the injected spin polarization 
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and the applied magnetic field lie along the same axis (y- 
axis) and the spin relaxation rate is given by longitudinal 
spin relaxation equation IV. 36 of Fabian et al. 12 The 
equation is rewritten using uj a = g e [iBAB exa /h and uiq = 



g e IJ-BBtotal/fi=9e^B(Ba 



B ex ,y/K) to yield, 



where the g* is the enhanced g-factor due to the presence 
of the exchange field. Substituting this into equation 
(6) yields the expression for spin relaxation from para- 
magnetic moments in the linear regime, 

1 (AS) 2 1 



1 



(ABl 



1 



(B, 



app,y 



B e 



,,) 2 +( 



(6) 



B„ 



(AB) 2 



(B a pp,y) 2 



where (AB) 2 = (AB ex ) 2 z + (AB ex ) 2 z . In other words, 
precession around randomly fluctuating exchange fields 
along the x- and z-axes induce spin relaxation. Equa- 
tion (6) also shows that the spin relaxation is suppressed 
by a large applied magnetic field. Intuitively, this occurs 
because the precession axis is defined by the large ap- 
plied field (along y-axis) and fluctuating fields along the 
x- and z-axes have very little ability to tilt the precession 
axis. This peak in spin relaxation at low magnetic fields 
produces the observed dip in Rnl- 

The presence of the average exchange field B ex ^ y in 
equation (6) shows that the spin relaxation is maximized w = gefJ-BB to tai/fi = 9e^B(B, 



(11) 



Thus, the longitudinal spin relaxation rate is a Lorentzian 
with a peak at zero applied field. 

For the Hanle geometry, the injected spin polar- 
ization is along the y-axis and the applied magnetic 
field is along the z-axis. In this case the spin re- 
laxation rate is given by transverse spin relaxation 
equation IV. 38 in Fabian et al. 12 The equation is 
rewritten using u a — g e fiBAB exa /h (for a = x,y,z) and 



app,z 



B ex , z )/h to yield 



when the B apPtV +B ex . y -- 
B P 



= 0, or B n 



-B exy . Because 



is proportional to the magnetization, paramag- 
netic moments will generate a dip in Rnl centered at 
zero applied field while ferromagnetic ordering will gen- 
erate a hysteretic dip centered away from zero applied 
field. The observed dip (Fig. le and If of main text) is 
centered at zero applied field and is not hysteretic, signi- 
fying the observed magnetic moments are paramagnetic. 

For paramagnetic moments, B exy takes the form of 
the Brillouin function (Bj) and is given by 



B, 



BM) 



VMA ex (SM,y) /geHB 

VMA ex JBj(^)/g e ii B 



2. J 



cothf — 

2J V 2J 



coth 

2J 



23 



(7) 



q (8) 



where J is the total angular momentum quantum 
number of the magnetic moment, £ = J k"^T Bapp.y, 
kg is Boltzmanns constant, and T is temperature. 
For our experiments at T=15 K, this reduces to 
£ = Jg e B app ,y/(22.32 Tesla). Thus, for the values of 
B a pp,y in our experiments £ -C 1 so that £?j«(J + 1)£/3 J 



to yield 



Bex,y — 



T) M A ex J(J + 1) 



B, 



app,y 



3fI B V 22 - 32 Tesla / 

Thus, the total field can be written as 



(9) 



-E$total — ^app.y H~ $ex,y 



-ex 
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(ABf 
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(12) 



where the fluctuating field is assumed to be isotropic: 
(AB ex ) 2 x = (AB ex ) 2 y = (AB ex ) 2 . For paramagnetic mo- 
ments, this becomes 



T 2 



(AB)- lg. 



M-Y 



(AB) 



(B a 



(13) 



IV. APPLICATION OF THE EXCHANGE 
FIELD MODEL TO EXPERIMENTAL R NL AND 
HANLE DATA 

The spin relaxation, diffusion coefficient, and interfa- 
cial spin polarization of the pristine sample are deter- 
mined though analysis of Rnl an d spin precession mea- 
surements. Fitting of spin precession data to the Hanle 
equation 



Rnl — S 



^-L 2 /4Dt 

V4ttD< 



,'ge^B D 

cos — — B n 



t)e 



-t/T> 



dt (14) 
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rj M A ex J(J + 1) 
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rj M A ex J(J + 1) 
3/x B (22.32 Tesla) 
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~B„ 



(10) 



provides values of spin lifetime (t so ), diffusion coeffi- 
cient (D), Hanle amplitude (5), and spin diffusion length 
(X = y/D t so ). For pristine graphene, the electron g- 
factor, q e , is assumed to be 2. It should be noted that 



G 



in the case of pristine graphene, the longitudinal spin re- 
laxation (rf°) and transverse spin relaxation (t|°) due 
to spin orbit coupling arc equivalent (Tf° = r|° = T ;io ) 14 . 
Data measured on sample A (Fig. Id of main text), yields 
t so —A79 ps, D=0.023 m 2 /s and A=3.3 pm for the channel 
length, L—5.25 /im. The corresponding Rnl data is fit 
with the non-local resistance equation (15) to obtain the 
interfacial spin polarization, Pj, for the graphene device. 




In the above equation, Rc is the graphene spin resis- 
tance defined by, Rg= Pg^/w, where pc is the resistivity 
and w is the graphene width. R\,2 denotes the contact 
resistances of injector and detector electrodes, Pf is the 
ferromagnetic electrode spin polarization (assumed to be 
0.35 for cobalt), and Rp = ppj^ is the spin resistance of 
the ferromagnet, where pf = 5.8x 10~ 8 fhn is the resistiv- 
ity of cobalt, lj=50 nm is the effective spin injector con- 
tact length of the ferromagnetic electrode and is deter- 
mined by the fabrication procedures (see Han et al. 15 for 
details) , and lastly, Xp is the spin diffusion length of the 
ferromagnet, taken to be 38 nm in cobalt. The measured 
Ai?ATL=8.8 57 for Sample A (Fig. lc) corresponds to a 
P/=0.20 for parameters A=3.3 /im, p G =898 CI, w=2.3 
pm, i?i=15.76 kfi, and i? 2 =4.00 kft. The contact re- 
sistances are measured in a three terminal geometry 16 
and are found to be unaffected by hydrogen exposure. 
The measured value of 20% interfacial spin polarization 
is comparable to previously reported values for efficient 
spin injection into SLG through tunneling contacts 16 . Pj 
is assumed to remain constant throughout hydrogen ex- 
posure, a reasonable assumption since the graphene at 
the site of spin injection is protected by the electrode, 
the hydrogen does not alter the cobalt (see section II C), 
and contact resistances remain unchanged. 

As discussed in the main text, exposure to atomic hy- 
drogen results in the formation of magnetic moments, de- 
tected as a dip in Rnl- Additionally, a sharpened Hanle 
curve signifies enhanced precession of injected spins due 
to the presence of an exchange field caused by the mo- 
ments. The exchange field is not accounted for in the 
standard Hanle equation (14), preventing direct deter- 
mination of t so and D through Hanle fitting. Instead, 
the Einstein relation 



D 



(16) 



is employed to obtain D for hydrogen-doped samples, 
where v denotes the density of states and e is the elec- 
tron charge. Assuming v is unchanged by exposure to hy- 
drogen, which is reasonable in the dilute limit, the diffu- 
sion coefficient of hydrogen-doped samples Dh y d is deter- 
mined from the pristine diffusion coefficient D pr i St i ne and 



the conductivities of hydrogen-doped (<Jhyd) and pristine 
{o-pristine) graphene. 



D 



hyd 



D 



pristine 



^pristine 



Vhyd 



6 V ' ^pristine 



(17) 



The change in conductivity from 1.113 mS to 0.143 
mS following hydrogen exposure results in Z)/ l j /c ;=0.0029 
m 2 /s. 

The longitudinal spin lifetime is evaluated by exam- 
ining Rnl for the hydrogen-doped sample. As shown 
in equation (11) the spin relaxation rate arising from 
the presence of magnetic moments is described by a 
Lorcntzian centered at B app y —Q, and can be fit using 
the general form 



1 



= r 



7 

(B ap p,y) 2 + 7 2 



(18) 



The total longitudinal spin lifetime, T{ otal , depends on 
both rf x and r so through the relation 



1 



I 1 



(19) 



subsequently affecting the spin diffusion length A 



A= \lDTi otal = \ l D ' ' ' 



<D 



(Bapp,y) 2 +J 2 



(20) 



The field dependent A directly translates to a field de- 
pendence in the non-local resistance causing the experi- 
mentally observed dip in Rnl at zero applied field. Val- 
ues for t so , r, and 7 are determined by fitting the mea- 
sured Rnl data to the non-local resistance equation (15), 
where A is field dependent and defined by equation (20). 
The best fit to sample A is obtained by r so =531 ps, 
r=2.73xl0 8 s"\ and 7=8.32 mT and is displayed as 
the red curve in Fig. 3a (main text), using, i?ir-=0.019 
ft, p G =6.99 kft i?i = 15.76 kO, i? 2 =4.00 kfl, Pj=0.20, 
Pp=0.35, and L=5.25 ^m. The field dependent values 
of Tf otal are displayed as the red curve in Fig. 3b (main 
text). Clearly, this model explains the data well and 
may also be relevant for dip features observed recently in 
metallic lateral spin valves 17 . 

Following the determination of r so , T, and 7, the spin 
precession data for the hydrogen-doped sample is ana- 
lyzed in order to obtain values for g*. The standard 
Hanle equation (14) must be modified to account for 
precession induced by both the applied field and the ex- 
change field [B e x,z] produced by magnetic moments as 
well as the field dependent transverse spin lifetime T\ otal , 
where Tt l tal = -b + i ■ For the Hanle geometry, the spin 



relaxation rate from the magnetic moments is given by 



7 



equation (13). Thus, in terms of the Lorentzian param- 
eters r and 7, the total spin relaxation rate is, 



1 



1 



rptotal 



7 



(B 



app , z ) 



■T 



(21) 



For comparison, Fig. 3b (main text) displays the to- 
tal longitudinal spin lifetime T* otal (red curve) and total 
transverse spin lifetime Tj ota/ (black curve) for hydrogen- 
doped graphene sample A. 

The modified Hanlc equation is dependent on T\ otal 
and takes the form 



Rnl — S 



e -L 2 /4Dt _ t 
y/A-KDt 6 



+ ir 1 + 





(9e^B lT > 

x cos(^ — — {B apPtZ 



B ex , z )t)dt 



(22) 



For paramagnetic moments, B ex z is described by the 
Brillouin function. Additionally, ( 2 2^2Tcsia) ^ ^ so * na * 
B ex ,z can be represented by the low field approximation 



Bex,z — 



rj M A ex J(J + 1) 
3^ B 



B, 



22.32 Tesla) 



app.z 



resulting in, 



Rnl — S 



x cos 







t \Bapp : 



r] M A ex J(J+l)/ B app , z 
3/xb V22.32 Tesla 



(23) 



(24) 



rff 



which simplifies to, 



(I 



y/4irDt 



with, 



.9e = Se 



T] M A ex J(J + 1) 



3^ B (22.32 Tesla) 



(25) 



(26) 



Fitting of precession data to equation (25) yields values 
for Hanle amplitude (S) and <?*. In the fit, fixed parame- 
ters are: r=2.73x 10 s s _1 and 7=8.32 mT, as determined 
by analysis of non-local resistance, -D=0.0029 m 2 /s from 
the Einstein relation, and L=5.25 /im. The best fit of 
Hanle data (presented in Fig. 3c (main text)) results 
in a g* value of 7.13. We take a reasonable value for 
the exchange coupling of A ex ~ 1 eV and the paramag- 
netic spin of J = 1/2 expected for the unpaired electrons 
due to hydrogen adatoms on graphene. Using the value 
g* =7. 13, we can independently estimate the fractional 
filling density of hydrogen induced magnetic moments to 
be r/M ~ 1 % using equation (26). This is in reasonable 



agreement with the order of magnitude estimate of 0.1% 
determined from the resistivity for 8 s hydrogen exposure 
to sample A. 

The correlation time (r c ) and the rms fluctuations in 
exchange field {AB) are determined by comparing equa- 
tions (11) and (18) to give 



1 



9*eHB 7 



(AB) 2 = 7 



(27) 



(28) 



Using 3* =7.13 from the Hanle fit and r=2.73xl0 8 s" 1 
and 7=8.32 mT from the non-local fit, we obtain val- 
ues of r c =192 ps and (AB) 2 =4.59x 1CT 5 T 2 (or AB = 
^{AB 



{AB e 



=6.78 mT). 



V. 



MAGNETIC MOMENTS GENERATED BY 
LATTICE VACANCIES 



We investigate the effect of lattice vacancy defects in 
graphene. Several theoretical works suggest the sim- 
ilarity of magnetism due to vacancies and hydrogen- 
doping 18 ' 19 , as both should create magnetic moments in 
graphene due to the removal/hybridization of p 2 -orbitals. 
It is therefore reasonable to expect that similar effects 
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FIG. 5. The effect of lattice vacancies on charge and spin 
transport in SLG at 15 K. a, Non-local spin transport data in 
pristine graphene with A_Ratl=38.9 Q. Hanle precession data 
(inset) yields values of r so = 859 ps and D=0.023 m 2 /s. b, 
Gate dependent pa for pristine graphene (black) and follow- 
ing lattice vacancy formation via Ar-sputtering (red), c, Non- 
local spin transport data in SLG containing vacancies. After 
Ar-sputtering, ARml is decreased to 8.2 Q and exhibits a dip 
at zero applied field. The minor loop (green) confirms the ob- 
served dip is a spin dependent effect, caused by paramagnetic 
moment formation. The best fit to the Hanle precession data 
(inset) yields g*=5.86. d, The field dependent T^otal used to 
fit the Hanle data in (c). 
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will be observable in graphene spin transport following 
the introduction of lattice vacancies. To induce vacan- 
cies on pristine SLG spin valves, in-situ Ar-sputtering is 
performed at a sample temperature of 15 K. Argon par- 
tial pressures of 1 x 10~ 6 torr and energies between 100 eV 
and 500 eV combined with short sputtering times (several 
seconds) produce dilute lattice vacancies. Prior to expo- 
sure to Ar-sputtering, the SLG device exhibits a ARnl 
of 38.9 Q at Vq — Vd=20 V and displays no dip in non- 
local resistance at zero applied field (Fig. 5a). Fitting 
of the corresponding precession data (inset of Fig. 5a) 
results in values of t so =859 ps and D=0.023 m 2 /s for the 
pristine SLG device. The black (red) curve presented in 
Fig. 5b displays pa before (after) sputtering. After the 
introduction of vacancies, the resistivity is substantially 
increased and the mobility is reduced from 4945 cm 2 /Vs 
to 949 cm 2 /Vs. Ar-sputtering results in a large decrease 
in the magnitude of ARnl as well as the emergence of 
a dip in Rml at zero applied field (Fig. 5c). The minor 
loop, shown in green, indicates the observed dip is due to 
a decrease in the spin signal at low fields, signifying the 
formation of paramagnetic moments. The Hanle data 
(Fig. 5c inset) narrows following Ar-sputtering. The 
Hanle data combined with fitting the dip in Rml yields 
values of 5* =5.86, AB=13.9 mT, r c =64.1 ps, and the 
field dependent T\ otal shown in Fig. 5d. 



VI. PROPERTIES OF THE EXCHANGE FIELD 

In this section, we discuss properties of the exchange 
field stemming from the formation of magnetic moments 
with the introduction of hydrogen. Specifically, we ex- 
amine the relation between the exchange field and nar- 
rowing of the Hanle curve, and we investigate the gate 
dependence and accuracy of g* . 



A. Exchange held and narrowing of the Hanle 
curve 

Conventional Hanle analysis, as described at the begin- 
ning of supplemental section IV, consists of fitting spin 
precession data to the Hanle equation (14), yielding val- 
ues for the spin lifetime (r so ), the diffusion coefficient 
(D), and the amplitude (S). This relies on the assump- 
tion, gl—2, and the absence of an exchange field. In 
conventional Hanle analysis, a narrowing of the Hanle 
curve is typically associated with an increase of the spin 
lifetime. Therefore, a valid question is whether the ob- 
served narrowing in the spin precession data after hydro- 
gen doping is due to an enhanced spin lifetime instead 
of the emergence of an exchange field. Comprehensive 
analysis comprising the full data set (conductivity, non- 
local spin resistance {Rnl), and Hanle spin precession) 
provides compelling evidence that the narrowing of the 
Hanle curve is due primarily to an exchange field as op- 
posed to enhanced spin lifetime. First, an increase in the 
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FIG. 6. Conductivity and calculated diffusion coefficient for 
pristine and hydrogen-doped SLG at 15 K. a, Gate depen- 
dent conductivity for pristine graphene (black curve) plotted 
on the left axis and the calculated diffusion coefficient, Dc 
(red curve), plotted on the right axis, b, Gate dependent 
conductivity for 8 s hydrogen exposure of SLG (black curve) 
plotted on the left axis and Dc (red curve) plotted on the 
right axis. 



spin lifetime cannot explain the observed dip in Rnl, 
while a fluctuating exchange field explains the dip and 
lineshape very well (see supplemental sections III and 
IV). Second, conventional Hanle analysis (g'*—2) of the 
hydrogen-doped sample A yields values for the diffusion 
coefficient that are inconsistent with the values obtained 
from the conductivity (differ by a factor of ~6) and val- 
ues for spin lifetime T\ otal that are inconsistent with the 
values obtained from the non-local spin resistance (differ 
by factor of 5-60). These inconsistencies can be resolved 
if an exchange field is present (g* > 2). Key features of 
the full data set emerge only after hydrogen doping and 
are best explained with a single effect, providing strong 
evidence for the presence of exchange fields. 

In the following, we provide a detailed analysis of the 
discussion outlined above. First, we investigate the Hanle 
spin lifetime, Tj ' ' —t so , without any consideration of 
an exchange field (i.e. g*=2). Fig 7a and 7c show the gate 
dependence of Hanle lifetimes obtained from fitting spin 
precession data using Hanle equation (14) for pristine 
and 8 s hydrogen exposure to sample A, respectively. As 
can be seen in Fig. 7c, when the Hanle fit parameters D 
and T2 are allowed to vary, best-fit values yield long spin 
lifetimes. Values of D from the Hanle fit are denoted as 
Ds and are displayed (black open squares) in Fig. 7b 
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FIG. 7. Conventional Hanle fitting assuming <?*=2 (no ex- 
change field), a, Gate dependence of the Hanle spin lifetime 
for pristine SLG obtained by the Hanle equation, b, Gate de- 
pendence of the diffusion coefficients Ds (black open squares), 
obtained from the Hanle equation and the calculated diffu- 
sion coefficient, Dc (red curve), obtained by calculation from 
the conductivity and the broadened DOS for pristine SLG. 
c, Gate dependence of the Hanle spin lifetime for 8 s hydro- 
gen exposure obtained by the conventional Hanle method, d, 
Gate dependence of the diffusion coefficients Ds (black open 
squares), obtained from the Hanle equation, and Dc (red 
curve) for 8 s hydrogen-doped SLG. e, Relative difference be- 
tween Ds and Dc for pristine (purple solid squares) and 8 s 
hydrogen-doped (light blue solid squares) graphene. 



and 7d for pristine and 8 s hydrogen-doped SLG. 

Alternatively, one can use the gate-dependent conduc- 
tivity to determine D via the Einstein relation (16), 

i> = 4- 

where a is the experimentally measured conductivity, e 
is the electron charge, and v is the broadened density of 
states (DOS). This value of diffusion coefficient is denoted 
as Dc- Spatial fluctuations of the Fermi level due to inho- 
mogeneities in the SiC^/Si substrate lead to broadening 
of the DOS. For Gaussian broadening 20 , the DOS is 



1/(25) 



h 2 v 2 F 



2b 



Serf 



E 

&7I 



(29) 



where g v is the valley g- factor, g s is the electron spin g- 
factor, h is Plancks constant, Ujr=lx 10 6 m/s is the Fermi 



velocity, and b is the Gaussian broadening parameter. In 
Figure 6a and 6b we show the conductivity (black curve) 
for both pristine and 8 s hydrogen exposure, respectively. 
These conductivity curves correspond with the resistance 
data of sample A as shown in Fig. lb (main text). We 
find reasonable agreement for 6=100 meV and use this 
throughout the remainder of this section for Dc- On the 
right axis of Figures 6a and 6b we plot the calculated 
diffusion coefficient (red curve) (&=100 meV) for pristine 
and 8 s hydrogen, respectively. We have found that a 
broadening parameter between 75 and 125 meV, which is 
reasonable for graphene on Si02/Si substrate 20 ' 21 , gives 
generally similar results for the present discussion. 

We next examine the difference between these two 
methods for determining the diffusion coefficient. Figure 
7b and 7d plots Ds (black open squares) and Dc (red 
curve) for the pristine sample and the hydrogen-doped 
sample as a function of gate voltage, respectively. Inter- 
estingly, Ds is much smaller than Dc for the hydrogen- 
doped sample, particularly when compared to the pris- 
tine sample. To quantify this, we plot the relative dif- 
ference (Dc — Ds)/Ds in Figure 7e and find it to be as 
large as ~6 for the hydrogen-doped sample. On the 
other hand, the relative difference is less than ~1 for the 
pristine sample. Therefore, for the hydrogen-doped sam- 
ple, the values of D determined from the conventional 
Hanle method (Ds) and the charge transport measure- 
ment (Dc) are inconsistent. There are two possible ex- 
planations for the appearance of a large discrepancy in 
Ds and Dc upon hydrogen doping. First, a system with 
an exchange field and increased effective g-factor yields a 
very low value of Ds when fit using conventional Hanle 
with 5g=2. Alternatively, it is well known that Ds and 
Dc can differ drastically if there are significant electron- 
electron interactions present in the system 22 . As dis- 
cussed below, we find that the presence of an exchange 
field also resolves other inconsistencies generated by con- 
ventional Hanle fitting. 

We now consider values of spin lifetime determined by 
the in-plane Rnl data and compare it to values deter- 
mined from conventional Hanle fitting assuming no ex- 
change field. Following the same procedure outlined in 
supplemental section IV from equation (15) to (20), we 
obtain values of r so based on best fits to the high field 
data of Rnl- The method utilized in section IV takes the 
diffusion coefficient, Ds, from the conventional Hanle fit- 
ting of the pristine sample, then scales it based on the 
Einstein relation according to equation (17). In this sec- 
tion, we denote this value as Dss- The resulting spin 
lifetime values from Rnl are plotted in Figure 7b as a 
function of gate voltage and labeled "£>ss". Alterna- 
tively, best-fit values for the gate dependence of spin life- 
time using the calculated diffusion coefficient, Dc, given 
by equations (16) and (29), are plotted in Figure 8b 
(blue open diamonds) and labeled "Dc" ■ The light blue 
shaded region in the range of 300-600 picoseconds repre- 
sents the values of spin lifetime consistent with the non- 
local resistance data and is labeled "Rnl" ■ We compare 
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FIG. 8. Comparison of the spin lifetimes, a, Gate depen- 
dence of the diffusion coefficients Ds (black open squares), 
Dc (red curve), and Dss (blue solid circles) for 8 s hydrogen- 
doped graphene. b, Gate dependence of spin lifetimes for 
graphene after hydrogen exposure. Grey shaded region rep- 
resents spin lifetimes by conventional Hanle fitting with Dc 
(black open diamonds) and Ds (black open squares) with- 
out taking the possibility of magnetic moments into account 
(<7e=2). A maximum T2 of 33 ns is obtained at the charge neu- 
trality point (Va= -1 V) for Dc=0.0032 m 2 /s. Blue shaded 
region represents spin lifetimes consistent with the in- plane 
Rnl data. Best fits solutions for t s ° using Dc (blue open 
diamonds) and Dss (blue solid circles) according to the pro- 
cedure outlined in supplemental section IV. 



this with the spin lifetimes determined by conventional 
Hanle fitting. The spin lifetime from the conventional 
Hanle fitting for 8 s hydrogen exposure displayed in Fig- 
ure 7c is replotted in Figure 8b and labeled "Ds" (black 
open squares). Alternative values for spin lifetime are 
obtained by performing the Hanle fit with the diffusion 
coefficient as a fixed parameter given by Dq (and g*=2). 
The resulting spin lifetime as a function of gate voltage 
is plotted in Figure 8b and labeled u Dc" (black open 
diamonds). The grey shaded region between 3 and 33 
nanoseconds represents the values of spin lifetime consis- 
tent with the Hanle data assuming g* —2 and is labeled 
11 Conventional Hanle" . Based on Figure 8b, the spin life- 
time determined by conventional Hanle analysis (with 
g* e =2) is inconsistent with the spin lifetimes determined 
by non-local resistance. 

To summarize, with the introduction of atomic hydro- 
gen (or lattice vacancies) to SLG, conventional Hanle fit- 
ting with the assumption of g* e =2 yields two inconsis- 
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FIG. 9. Gate dependence of the enhanced precession due 
to the exchange field, a, Hanle precession data at the charge 
neutrality point (Vg=-1 V) for 8 s hydrogen exposure at T=15 
K. A sharp increase in the observed precession due to the 
exchange field is seen at the charge neutrality point. Best 
fit to the Hanle equation (red curve) taking into account the 
Lorentzian shape of the spin lifetime and the scaled diffusion 
coefficient, Dss, yields an effective g-factor, p*=19.2. b, Gate 
dependence of g* for both Dc (blue open diamonds) and Dss 
(blue solid circles). 

tencies: (i) values of Ds that are improbably low when 
compared to Dc and (ii) values of spin lifetime that are 
too large compared to values obtained from the non-local 
resistance. Notably, both of these inconsistencies can be 
alleviated if g* > 2. This can be understood by consider- 
ing the symmetries of the Hanle equation (30), 

OO r. 

/p—li/iDst /n*u \ 

—— Cos i9eM Ba A e -t/T 2dt (30) 



This equation is invariant under the transformation 
g* — ► eg*, T-2 — > T2/C, Ds — > cDs, S — > cS, where c is a 
constant. For a given parameter set (g*, T2, Ds, S), the 
transformed parameter set (cg*,T 2 /c,cDs,c S) will gen- 
erate the same Hanle curve. Therefore, if we begin with a 
conventional Hanle fit that assumes g* e =2 1 a transforma- 
tion that increases g* (i.e. c> 1) has the effect of decreas- 
ing T 2 and increasing Ds- This simultaneously alleviates 
the discrepancies in spin lifetime and diffusion coefficient 
mentioned earlier, and therefore provides strong evidence 
for the presence of exchange fields. 

The discussion above highlights several key points 
about the analysis of Hanle data. In the presence of 
an exchange field, g* e becomes a free parameter and the 
fitting parameters cannot be determined uniquely from 
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the Hanle data alone. Therefore, it becomes necessary 
to analyze the in-plane Rnl data and the Hanle data 
together (as detailed in supplemental section IV) in or- 
der to determine key parameters such as spin lifetime 
and g*. The nature of this data set, with in-situ dop- 
ing, makes it straightforward to apply this analysis, but 
this may not be true in other studies utilizing Hanle spin 
precession. Consequently, it brings to light an important 
question about the use of Hanle fitting in general: how 
does one tell whether changes in the Hanle curve are due 
to changes in spin lifetime or g*7 Fortunately, the above 
analysis leads to a useful rule of thumb: If values of Dc 
and Dg are similar for conventional Hanle fitting, this 
provides support for the assumption that g%=2. This is 
important for future Hanle studies of spin relaxation in 
order to recognize when changes in spin precession data 
are due to changes in spin lifetime. For systems with an 
exchange field, this analysis motivates the need for alter- 
native experimental techniques that can independently 
measure g* e and T2, such as electrically-detected electron 
spin resonance (ESR) and time resolved spectroscopies. 



exchange field. Following the procedure of supplemental 
section IV, g* values are obtained by Hanle fits to the 
spin precession data for 8 s hydrogen exposure to sample 
A using equation 25. Figure 9a shows the spin preces- 
sion data (black closed circles) at the charge neutrality 
point (CNP) (Vg=-1 V) where the fastest precession is 
observed. The best fit solution (Fig. 9a red curve) to 
the spin precession data by equation 25 is determined 
through the free parameters S and g*. Fig. 9b displays 
the gate dependence of g\ for both D$s (blue solid cir- 
cles) and Dc (blue open diamonds) . The minimum value 
is g*=5.1 obtained at Vq—-1\ V and the maximum is 
<7*=19.2 at the CNP assuming D ss . For D Cl the mini- 
mum g* value is 10.3 for Vg= -21 V and the maximum is 
<7*=21.1 at the CNP. Uncertainty of the D value leads to 
uncertainty in g*, again highlighting the need for tech- 
niques in graphene spintronics to directly measure g* and 
T 2 . 



B. Gate dependence and accuracy of g% 

In this section we examine the gate dependence of the 
effective electron g-factor, g*, due to the presence of an 
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